Draft version February 5, 2008 

Preprint typeset using I£T^X style cmulateapj v. 10/09/06 



O 
O 

(N 

(N 
> 

On 

t> 

O 

Oh' 
i 

O 

CO 

c3 



X 



CONSTRAINTS ON NEUTRON STAR PROPERTIES FROM X-RAY OBSERVATIONS 

OF MILLISECOND PULSARS 

Slavko Bogdanov, George B. Rybicki, and Jonathan E. Grindlay 

sbogdanov@cfa.harvard.edu, grybicki@cfa.harvard.edu, jgrindlay@cfa.harvard.edu 
Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 
Draft version February 5, 2008 

ABSTRACT 

We present a model of thermal X-ray emission from hot spots on the surface of a rotating compact 
star with an unmagnetized light-element atmosphere. An application to ROSAT, Chandra, and XMM- 
Newton X-ray observations of the nearest known rotation-powered millisecond pulsar (MSP) PSR 
J0437-4715 reveals that the thermal emission from this pulsar is fully consistent with such a model, 
enabling constraints on important properties of the underlying neutron star. We confirm that the 
observed thermal X-ray pulsations from J0437-4715 are incompatible with blackbody emission and 
require the presence of an optically thick, light element (most likely hydrogen) atmosphere on the 
neutron star surface. The morphology of the X-ray pulse profile is consistent with a global dipolc 
configuration of the pulsar magnetic field but suggests an off-center magnetic axis, with a displacement 
of 0.8 — 3 km from the stellar center. For an assumed mass of 1.4 M Q , the model restricts the allowed 
stellar radii to R = 6.8 — 13.8 km (90% confidence) and R > 6.7 km (99.9% confidence), which is 
consistent with standard NS equations of state and rules out an ultracompact star smaller than its 
photon sphere. Deeper spectroscopic and timing observations of this and other nearby radio MSPs 
with current and future X-ray facilities {Constellation- X and XEUS) can provide further insight into 
the fundamental properties of neutron stars. 

Subject headings: pulsars: general — pulsars: individual (PSR J0437-4715) — stars: neutron 
X-rays: stars — gravitation — relativity 



1. INTRODUCTION 

After four decades of considerable observational and 
theoretical investigation of neutron stars (NSs), very lit- 
tle is known about several key properties of these objects. 
In particular, there is very limited information concern- 
ing their magnetic field configuration, surface properties, 
and structure and composition of their interiors. Ther- 
mal radiation from the physical surface of a NS can po- 
tentially serve as a useful tool in the study of these prop- 
erties. For instance, a measure of the surface tempera- 
ture can provide insight into the structure of the stellar 
interior through comparison with theoretical NS cooling 
models (see, e.g., Page 1998, and references therein), al- 
though it relics on an accurate age determination which 
is not easily measurable for most systems. In addition, 
mapping the temperature variations accross the surface 
may be used to deduce the topology of the surface mag- 
netic field of strongly magnetized NSs. Surface emission 
can also be used to infer the compactness of the star 
through measurement of the gravitational redshift z g . 
In practice, this endeavor has proven to be quite diffi- 
cult due to the abs ence of spectral lines (see, however, 
ICottam et al.l 12006). For a NS radiating from the entire 
surface, the gravitational redshift can be obtained by de- 
termining the effective emission radius, i? e ff- However, 
this method assumes that the inferred radius is exactly 
equal to the radius at infinity, R°° — (1 + z g )R, of the 
star. As such it is subject to uncertainties in the the 
emission properties of and exact temperature distribu- 
tion across the stellar surface, as well as the distance to 
the source. 

A promising approach is to study systems in which the 
thermal radiation is confined to a very small fraction of 



the stellar surface (R e g <C R). Such emission geometry is 
observed in rotation- and accretion-powered milisecond 
pulsars (MSPs) as well as old normal pulsars. In these 
objects, modeling the spectrum and rotation-induced 
variations of the observed flux c an provide a measure 
of the compactness of the NS (see lPavlov fe ZavTm1ll997t 
IZavlin fc Pavlovi[l998h . Moreover, the hot spots on the 
surface are believed to coincide with the polar caps, 
where the magnetic field is anchored to the NS surface, 
thus, allowing one to determine the geometry of the NS 
magnetic field. In addition, the modulation of the ob- 
served flux caused by the changing projection of the emis- 
sion area permits study of the radiation pattern of the 
emergent intensity (isotropic versus anisotropic). 

Herein, we present model spectra and lightcurves of 
compact stars with hot spots covered by a hydrogen 
atmosphere. An application of our model to archival 
ROSAT, Chandra, and XMM-Newton observations of 
the radio MSP J0437-0437 permits us to gain valu- 
able insight into several fundamental NS parameters. 
Th is work represents an exte nsion of similar studie s 
bv lPavlov fc Zavlinl (|1997f ) and IZavlin fc Pavlov! i|1998t ). 
The present paper is organized in the following manner. 
In §2 we outline the theoretical model, while in §3 we 
describe its practical application. We offer a discussion 
in §4 and give conclusions in §5. 



2. THEORETICAL MODEL 
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120041 : ICadeau et alj|2007t) . However, with some notable 



exceptions (e.g., IZavlin et al.l Il995at IPavlov fe ZavlinI 



Il997t IZane fc Turollall2006l ). these studies have assumed 
the idealized case of blackbody emission, i.e. an emer- 
gent intensity independent of angle. Here, we examine 
in-depth the observable properties of compact stars 
with light-element atmospheres and illustrate the key 
differences from blackbody emission. Below, we outline 
the basic formalism describing the system geometry and 
photon propagation. 

2.1. Geometry 

The basis of the model consists of a compact star of 
mass M, radius R, spin period P, and two antipodal 
hot spots 1 with i? c ff <C R. We designate the hot spot 
that approaches closer to the observer as the primary 
and define the closest approach as rotational phase zero 
(</> = 0). The magnetic axis is inclined at an angle a 
with respect to the spin axis, while the spin axis is at an 
angle £ relative to the line of sight to a distant observer. 
The position of the primary hot spot on the NS surface 
is defined by the angle between the normal to the surface 
and the line of sight: 

cosip(t) = sin a sin £ cos <f>(t) + cos a cos £ (1) 

An expression for the secondary hot spot is easily ob- 
tained by substituting a — > ir — a and <f> — * % + 4>. Note 
that a and £ are interchangeable as they figure in equa- 
tion (1) in the same way. 

In Schwarzschild geometry, due to the compact na- 
ture of the NS, the gravitational field has a profound 
effect on the photons as they propagate from the stellar 
surface to infinity. In particular, the energy of a pho- 
ton is reduced by a factor of 1 + z g = (1 — Rs/R)^ 1 ^ 2 , 
where Rs = 2GM/c 2 7 as it escapes from the deep gravi- 
tatio nal potential (see, e.g.. iMisner. Thorne. fc Wheeler! 
1970). Furthermore, a photon emitted at an angle 
9 > with respect to the local radial direction follows a 
curved trajectory and is observed at infinity at an angle 
i/j > 9. The relation betw een these two angles is given 
by (|Pechenick et al.lll983f) : 



where 



dr 
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1 
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Rs 



-1/2 



y/1 - Rs/R 



(2) 



(3) 



is the impact parameter at infinity of a photon emitted 
from radius R at an angle 9. In many practical situa- 
tions, one can use a m ore convenient approximate rela- 
tion between ip and 9 (jZavlin et al.lll995at iBeloborodovl 
[2002T) 

cos9-R s /R 
C ° S ^ I -Rs/R (4) 

which is valid for R > 2R$ and is remarkably accurate 
(fractional error of order a few percent). One interesting 
consequence of the bent photon trajectories is that most 
of the NS surface is observable at any given time. In 

1 In §3 we also consider the case when the hot spots are not 
antipodal. 



flat spacetime, the visibility condition is simply cos?/' = 
cos 9 > 0, while in Schwarzschild geometry a point on 
the neutron star surface is visible down to a critical angle 
cosi/> c , corresponding to fe max = R°° = R/ y/1 — Rs/R, 
the apparent radius of the star. 

2.2. Doppler Effect and Time Delays 

For rotation- and accretion powered MSPs, the rapid 
motion of the NS surface induces an appreciable Doppler 
effect, parameterized through the familiar Doppler factor 

V = — (5) 



7(1 — v/ccost;) 



where 



7 



1/V1 - {v/cf, V 



2wR/P(l 



i?s/i?) -1 / 2 sin a is the velocity of the hot spot as mea- 
sured in the inertial frame of the hot spot, and £ is the 
angle between the direction of t he velocity vector and 
the lin e of sight. As shown by [Poutanen fe Gierlinskil 
(|2003h and lViironen fc PoutanerJ (|2004h . this angle can 
be expressed in terms of 9, ip, o>, and <fi as 



sm9 

cos 4 — : — - sin a sin < 

smip 



(6) 



where the quantity sin 9/ sin-0 can be approximated by 
its asymptotic value for small angles (1 — Rs/R) 1 / 2 - 

We note that the assumption of a spherical Kerr metric 
does not lead to an appreciable difference in the shape of 
the lightcurves c ompared to the Sch warzschild case (as 
demonstrated by ICadeau et al.l [20071. We also empha- 
size that for P < 3 ms deviations from spherical sym- 
metry su ch as rotation- i nduce d oblateness becomes im- 
portant. Ca deau et al.l (|2007f ) have found that in that 
regime neither the spherical Schwarzschild nor spherical 
Kerr metrics provide an accurate description of the prop- 
erties of the rotating star. Therefore, the validity of our 
model is limited to spin periods greater than ^3 ms. 

Photons emitted from the far side of the NS relative to 
the observer, in addition to following a curved trajectory, 
have to travel a greater distance compared to a photon 
emitted radially. The resulting time lag of the photon 
as seen by a distant observer is given by the elliptical 
integral (Pcchcnic k et al.|[l983h 



At(b) 



dr j 

1-R S /R\ 



Rs 



-1/2 



(7) 

This time delay translates into a phase lag (A</>) of a 
photon 

2?r A 
Acf> = —At 



P 



(8) 



which yields the observed phase </> b s = (/) + Acf) 
(|Viironen fc Poutanenl l2004). For R/R s « 2.5, the max- 
imum At, corresponding to a photon with maximum im- 
pact parameter 6 max = (1 + z g )R, is of order 60 /is, or 
<6% of the pulsar period for P > 1 ms. Therefore, 
propagation time differences have a minor effect on the 
observed lightcurves even for the most rapidly rotating 
neutron stars. Nonetheless, this effect has been included 
in our model for completeness. 

The observed flux per unit frequency from each spot is 
given by 

F(v) = I(v)dn (9) 
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Fig. 1. — Representative model lightcurves for a rotating M = 1.4 
Mq, R = 10 km NS with two point-like antipodal hot spots. The 
panels show classes I (a = 10°, £ = 30°), II (a = 30°, C = 60°), 
III (a = 60°, C = 80°), and IV (a = 20°, C = 80°) as defined 
by Beloborodov (2002). The solid lines in each plot correspond 
to a hydrogen atmosphere (black) and blackbody emission (grey) 
with no Doppler effect included. The dashed lines show the effect 
of Doppler boosting and aberration for P = 4 ms. All fluxes are 
normalized to the value for a = £ = 0. Two rotational cycles are 
shown for clarity. 

where I(v) is the intensity of the radiation as measured 
at infinity and dQ is the apparent solid angle subtended 
by the hot spot. Transforming both quantities to the rest 
frame of the hot spot yields 

F{v) = (1- Rs/R) 1/2 ^I\v\0>) cos9' (10) 

dcosi/) D z 

where the primed quantities are measured in t he NS 
surface rest frame (jPoutanen fc Gierlihskil [2003) . with 
cos<9' = r] cos and dScosfl = dS"cos6>'. I'(y',9') is the 
emergent intensity, dS' is the emission area and D is the 
distance. The three Doppler factors arise from the trans- 
formation of the intensity. A fourth factor is obtained 
upon integration of equation (10) over a frequency inter- 
val since dv = (1 - R s /R) 1/2 r]dv'. The total observed 
flux for a given rotational phase is found by relating </> 
and 9 for each hot spot through ip via equations (1) and 
(2), using the corresponding I'(v',6') in equation (10), 
and summing the flux from the two hot spots. Equation 
(10) allows one to generate a phase-integrated spectrum 
by simply integrating over the whole rotation period. It 
can also be used to construct an arbitrary emission re- 
gion on the NS surface by considering multiple surface 
elements. 

2.3. Atmosphere Model 

To describe the emission properties of the stellar 
surface we use the hydrogen atmosphere model from 



iMcClintock et al l (|2004h (see also iHeinke et all l2006h 2 
having the following set of assumptions. The atmosphere 
is static, in radiative equilibrium, and composed purely of 
hydrogen. The scale heights in the atmosphere are much 
smaller than the NS radius so a plane-parallel approxima- 
tion is valid. The star is weakly magnetized (B -C 10 9 
G), meaning that the effects of the magnetic field on 
the opacity and equa tion of state of the atmosphere can 
be ign ored (see e.g., IShibanov et al.l fT992t IZavlin et alJ 
Il995bl for a treatment of magnetized atmospheres) . The 
opacity within the atmosphere is due to thermal free- 
free absorption plus Thomson scattering in the unpo- 
larized, isotropic approximation. For the temperatures 
of interest (~10 6 K) we expect complete ionization im- 
plying that bound-free and bound-bound transitions are 
unimportant. Finally, this model is valid for tempera- 
tures below T e g w 3 x 10 6 K, where the effects of Comp- 
tonization within the atmosphere are neglegible. 

An important characteristic of light-clement NS atmo- 
sphere models is that the peak emission occurs at higher 
energies than a black body (BB hereafter) for the sam e 
effective temperature (|Romanilll987t [Zavlin et al.lll996f) . 
Furthermore, the emission pattern of the atmosphere 
is inherently anisotropic with intensity decreasing as a 
function of angle with respect to the normal, resulting i n 
a limb-darkening effect (cf Fig. 7 of IZavlin et al.l 11996). 
This implies that although the emission spectrum of such 
an atmosphere is qualitatively similar to the case of a BB, 
the observed rotation-induced variation of the hot spot 
flux will be profoundly different. Figure 1 shows a sam- 
ple of lightcurves, generated using equation (10) for a 
range of a and £ for both BB and H atmosphere models, 
assuming M — 1.4 M© and R = 10 km. The most no- 
table observable property of the atmosphere lightcurves 
is the substantially larger "depth" (i.e. pulsed fraction) 
relative to the BB case. Effectively, the limb darkening 
of the atmosphere acts to partially negate the effect of 
light bending, which tends to reduce the degree of modu- 
lation. Note, however, that even in the case of point-like 
hot spots, the modulations are still quite broad. In addi- 
tion, it is not possible to achieve arbitrarily large pulsed 
fractions (~100%) for any combination of a and £. These 
characteristics can be used to distinguish between ther- 
mal and non-thermal pulsed emission in cases where the 
spectral continuum alone does not provide sufficient in- 
formation (see e.g. IChatteriee et alJ 120071 for the case 
of PSR J0737-3039A). Limb darkening also significantly 
diminishes the effect of Doppler boosting on the shape of 
the pulse profile relative to the BB model, as evident in 
Figure 1. This is because, for a given hot spot, when the 
component of the velocity vector along the line of sight 
is largest (corresponding to a maximum Doppler boost) 
the hot spot provides a very small contribution to the 
total flux. 

Unlike a BB, the pulsed fraction of the atmosphere 
emission is expected to vary depending on the choice of 
energy band due to th e energy-dependen t optical depth 
of the hydrogen layer ([Zavlin et al.l [l996). This marked 
difference between the atmosphere and BB lightcurves al- 
lows one to discriminate between the two emission mod- 
els. For phase-resolved spectra of an atmosphere we ex- 

2 This model is virtually identical to that presented by 
IZavlin et all lll"996h . 
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Fig. 2. — Phase-integrated synthetic H atmosphere hot spot 
emission spectra for M = 1.4 M Q , R = 10 km, and T off = 2 X 10 6 
K. The solid lines correspond to classes I, III, II, and IV, from 
top to bottom, respectively, with values of a and £ as in Figure 
1. The dotted line is for a star with one hot spot always face-on 
(a = £ = 0). For all spectra the same arbitrary i? e g is assumed. 
Note the shifts in peak energy. 



pect softening of the emission at the pulse minima, due 
to the presence of a temperature gradient within the at- 
mosphere (with the coolest layer at the surface). For 
many systems, however, only phase-integrated spectra 
are available due to inadequate photon statistics and/or 
time resolution. Such model X-ray spectra in which the 
geometry, strong gravity, and rotation of the NS have 
been taken into account are shown in Figure 2. It is 
apparent that in addition to a change in the total flux, 
there is a significant shift in the peak energy of the spec- 
trum. For a BB spectrum no such variations are ex- 
pected. Therefore, even in the case of phase-integrated 
spectra it is important to take into account the system 
geometry (a and £) to ensure reliable measurements of 
the temperature and effective area of the hot spots. 

3. APPLICATION 

Close examination of equations (2) and (10) reveals 
that the flux observed at infinity from the hot spots de- 
pends on several important NS parameters. Figure 3 
illustrates the effect of the stellar compactness on the 
pulse profile of a rotating NS with two antipodal hot 
spots covered with a H atmosphere. It is apparent that 
a small increase in M/R results in a substantial decrease 
in the pulsed fraction. This is a consequence of the 
strong sensitivity of light bend ing on the comp a ctness 
of t he star. Thus, as show n bv I Pavlov fc Zavlml (|1997l ) 
and lZavlin fc Pavlov! (|1998ft . modeling of the spectra and 
pulse profile may, in principle, allow stringent constraints 
on M and R. Note that in the case of atmosphere models, 
M and R also figure implicitly in the intensity I'{v' , 9') 
through the surface gravity g s = (1 — Rs / 'i?) _1 / 2 A/ / R 2 . 
Although this dependence is much weaker than that 
through equation (2), it suggests that M and R should 
be treated as separate parameters. The spectra and 
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Fig. 3. — Model H atmosphere lightcurves for different stellar 
radii for an assumed 1.4 Mq star. The curves correspond to R = 9 
km (dot- dashed), 12 km (solid), and 16 km (dashed). The values 
of a and f for each panel are the same as in Fig. 1. Two rotational 
cycles are shown for clarity. 

lightcurves also depend strongly on the angles a and £ 
through equations (1) and (2). In rotation-powered pul- 
sars, a corresponds to the pulsar obliquity, i.e. the angle 
between the spin and magnetic axes. Knowledge of this 
angle has important implications in the study of NS mag- 
netic fields. Hence, an application of our model to X-ray 
spectroscopic and timing observations may yield useful 
constraints on M/R, a, and £. 

To allow direct comparison of our model with obser- 
vational data we have generated a grid of H atmosphere 
emergent intensities, I'(v',6'), for a range of effective 
temperatures T c g, g s (implying a range of M and R), 
and 6. This grid allows us to generate a spectrum and 
pulse profile for any desired values of M, R, a, and £. To 
minimize the number of free parameters, we fix M to 1.4 
M and vary only R. For the purposes of spectral analy- 
sis, we have incorporated our model int o XSPEC 3 12.3.0. 
The fi ts to the data are performed as in iPavlov fc Zavlinl 
(11997ft and lZavlin fc Pavlovl (|1998t) . We first fit the spec- 
trum for assumed values of R, a, and £. The resulting 
best fit T e ff, R e s (the effective radius of each hot spot), 
and Nh (the integrated H column density along the line 
of sight) are then used to generate the corresponding 
pulse profiles. After convolving with the appropriate 
detector response and accounting for background noise, 
the model lightcurves are binned to the same number of 
phase bins as the data, which allows a calculation of x 2 



X 



K 

E 



N h 



(11) 



where N Q ^, N m ^, and Nj,^ are the observed, model, and 

3 http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 
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TABLE 1 

X-ray Observations of PSR J0437-4715 



Observatory 


Detector 


Energy band 


Exposure 


Note a 






(keV) 


(fa) 




ROSAT 


PSPC 


0.1-2 


9.8 


S+T 


Chandra 


ACIS-S 


0.3-8 


25.7 


S 




HRC-S 


0.1-10 


19.6 


T 


XMM-Newton 


EPIC-MOS 


0.3-10 


68.3 


S 




EPIC-pn 


0.3-10 


67.2 


S+T 



References. — Zavlin & Pavlov 1998; Zavlin et al. 2002; 
Zavlin 2006. 

a S=spectroscopy, T=timirig. 



background counts in the ith phase bin, respectively, and 
K is the total number of bins. This procedure is repeated 
for a range of R, a, C, and JVh to map the multidimen- 
sional x 2 surface and derive confidence levels for these 
parameters. 

3.1. The Radio Millisecond Pulsar PSR 30437-4715 

Recent X-ray observations have revealed that 
the emission from most rotation-powered MSPs 
(iBecker k Ach cnbach 2002; iGrindlav et alll2002t IZavlinl 
l2006t Bogdan ov et al.l f2006af ) is of predominantly ther- 
mal nature. This radiation is believed to originate 
from the magnetic polar caps of the pulsar that are 
heated by a back-flow of high-e nergy particles from 
the magnetosphere (see e.g., lHarding k Muslimovl 
20020. Given that the progenitors of radio MSPs 
undergo an extended period of a ccretion duri n g the 
low-mass X-ray binary phase (|Alpar et al.l Il982t 
iBhattacharva k van den Heuvell I1991D it is reasonable 
to expect these stars to have accumulated a substantial 
surface layer of gas. Due to the immense surface 
gravity the accreted material quickly stratifies, with 
the lightest element remaining on top and dominating 
the observed emission. Therefore, the existence of a 
thick light elem ent atmosphere at t he surface of MSPs 
is highly likely (| Zavlin et al.l Il995af) . Note that for the 
quality of the available spectral data of MSPs, at the 
temperatures relevant for these sources a He atmosphere 
is virtually indistinguishable from a H one, although 
based on the standard formation scenario of MSPs the 
existence of H at the surfac e is m ore probable . Also , 
iRajagopal k Romanil (|1996l ) and IZavlin et all (|2002l ) 
have shown that a heavy element (Fe) atmosphere does 
not provide a good description of MSP spectra. 

Of the ~40 X-ray d etected radio MSPs, at present, 
only PSR J0437-4715 dJohnston et al.l ll993l. the near- 
est and brightest, has X-ray data of sufficient quality to 
per mit detailed modelin g of i ts spect r um a nd pulse pro- 
file. IZavlin et all (|2002h and IZavlinl (|2006h have shown 
that the spectrum of this pulsar is well described by a 
multi-component model consisting of two dominant hy- 
drogen atmosphere thermal components with tempera- 
tures Ti = (2.1-1.5)xl0 6 KandT 2 = (0.52-0.54) x 10 6 
K and effective radii i?i w 0.34 km and i?2 ~ 2.0 km, for 
assumed M = 1.4 M Q , R = 10 km, a = 45° and C = 45°. 
The presence of two temperatures can be attributed to 
non-uniform heating of the polar c aps by the return flow 
of magnetospheric particle s (see lHarding k Muslimovl 
12001 iZhang k Chendl2003h . 

A faint power-law (PL) tail is also observed from this 
MSP at energies >2 keV. The nature of this PL emis- 
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Fig. 4.— {Upper panels) ROSAT PSPC, Chandra ACIS-S, and 
XMM-Newton EPIC-MOS1/2 and EPIC-pn spectra of PSR J0437- 
4715. The lines show a representative best fit of our model to the 
data. The bottom panel shows the best fit residuals. (Lower panel) 
Representative best-fit model spectrum consisting of two hydrogen 
atmosphere and one PL components. 

sion is ambiguous as it may arise either due to mag- 
netospheric emission processes or weak Comptonization 
of the thermal radiation by the tenuous magnetospheric 
plasma (Bogdanov et al. 2006b). A crucial difference 
between the two interpretations is that in the former the 
PL component extends to energies well below the soft 
X-ray band, while in the latter the PL is only present at 
higher energies (>2 keV). Moreover, the magnetospheric 
emissi on PL is inconsistent w ith the FUV flux from this 
MSP (|Kargaltsev et al.ll2004l ) unless its spectral photon 
index is T < 1.6. 

Previosly, iPavlov k Zavlinl (|1997f ) and 

IZavlin k Pavlovl (|1998| ) have applied a two-temperature 
hydrogen atmosphere polar cap model to spectroscopic 
and timing ROSAT PSPC observations of PSR J0437- 
4715 in an effort to constrain M/R of this NS. In this 
analysis, a centered dipole field was considered. Based 
on this assumption a limit of R > 8.8 km for M = 1.4 
Mq was obtained. As deeper observations of this MSP 
have become available since, it is worth revisiting this 
object to determine if better constraints on M/R can be 
obtained. 

We have retrieved archival ROSAT, Chandra, and 
XMM-Newton spectral and timing observations of this 
MSP to investigate its X-ray emission. These obser- 
vations are summarized in Table 1. It is important to 
note that the values for the effective tem peratur es and 
radii derived by IZavlin et ail (|2002t ) and IZavlinl (|2006l ) 
are by no means unique as they are sensitive to the ini- 
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tial choice of M, R, and (as evident in Fig. 2) a and £. 
We have carried out a series of spectral fits for M = 1.4, 
R = 8 - 16 km, a = 0° - 90°, C = 0° - 90° 4 , and 
iVn = (1 — 3) x 10 19 cm -2 . For this set of parame- 
ters, we find T x = 1.43 - 1.85 MK, Ri = 0.04 - 0.40 km, 
T 2 = 0.4-0.56 MK, and R 2 = 1.8-4.9 km, with xl = 1-3 
for 402 degrees of freedom. The ranges quoted give the 
±1<T bounds for one interesting parameter. Given the 
known cross-calibration uncertainties of the various de- 
tectors, the minimum value of % 2 corresponds to an ac- 
ceptable fit and indicates a systematic uncertainty at the 
level of ~5-10%. These results are generally in good 
agreement with those derived in previous studies. We 
have performed the spectral fits for two cases: PL emis- 
sion with r < 1.6 and no PL emission below 2 keV as 
an approximation to a Comptonization model. In both 
cases, we find that such a component contributes with 
<5% to the total photon flux below 2 keV. Furthermore, 
although the statistics at high energies are quite lim- 
ited, XMM-Newton EPIC-pn data above 3 keV suggest 
th at the PL does not sho w sharp pulsations (cf Fig. 3 
of iBogdanov et alll2006bl ) as o bserved in the more lumi- 
nous MSP s PSRs B1937+21 (iNicastro et"aTl l2004D and 
B1821-24 (jRutledee etHl I2004D . This implies that the 
observed pulse shape is not significantly affected by the 
PL emission. Therefore, in our lightcurve analysis we will 
include this component as a constant (DC) contribution 
to the pulsed flux. 

3.2. Fits to the Pulse Profile of PSR 30^31-^115 

Using the derived range of temperatures and radii we 
proceed to generate synthetic pulse profiles and compare 
them to those of J0437-4715. When folded at the spin 
period (P = 5.76 ms), the X-ray emission from this MSP 
exhibits a single broad pulse, with pulsed fracti on 30- 
40%, depending on the choice of energy band fsee lZavlinl 
2006), consistent with thermal radiation. A visual com- 
parison with the lightcurves in Figure 1 shows that the 
X-ray pulsations of this MSP (Fig. 5) most closely resem- 
ble those of classes I and II. However, it is immediately 
apparent that the pulses of J0437-4715 are asymmetric, 
characterized by a faster decrease than increase in flux. 
This is clearly inconsistent with the antipodal hot spot 
model. In addition, at the spin period of this MSP the 
Doppler effect is fairly weak so it cannot reproduce the 
degree of asymmetry observed. Indeed, an acceptable fit 
could not be obtained for any combination of R, a, and 
(. The discrepancy can be resolved if the two hot spots 
are not diametrically opposite, i.e. if the magnetic axis 
is off-center. Therefore, in our model we introduce two 
additional free parameters, A<p and Aa, offsets in <j> and 
a, respectively, of the secondary hot spot. These corre- 
spond to displacements on the NS surface in longitude 
and latitude, respectively, from the antipodal position. 
The net offset on the stellar surface of the secondary hot 
spot from the antipodal position is 

As = i?cos~ 1 [— sinasin(a-l-Aa)— cos a cos(a+Aa) cosA</>] 

(12) 

4 For either a or f, the range 90° — 180° becomes equivalent to 
90° — 0° if we reassign the secondary hot spot to be the primary 



In an offset dipole configuration a is no longer the angle 
between the spin and magnetic axes and the dipole axis 
does not run through the center of the star. The total 
displacement (i.e. impact parameter) of the magnetic 
axis from the center is then 



,'As 

A , = Rsm I — 



(13) 



and apply the transformation 
necessary to consider the latter range. 



7r. Therefore, it is only 



If we introduce an offset dipole, we find that a phase lag 
of Acf) ~ —20°, implying Aa; ~ 1 km, is able to account 
for the asymmetry in the pulse profile. 

In the formal fits to the X-ray pulse profiles of PSR 
J0437-4715 we consider nine free parameters: T\, Ri, 
T2, i?2, -R, ot, £, Aa, and A<j>, while keeping JVh fixed 
at 2 x 10 19 cm -2 and M at 1.4 M . In addition, since 
the relative phase between the various observations can- 
not be determined due to the inadequate absolute timing 
precision of ROSAT and XMM-Newton, we allow the 
<p to vary independently for each dataset. The model 
was fitted individually to ROSAT PSPC (0.1-2.4 keV), 
Chandra HRC-S (0.1-10 keV), and XMM-Newton (0.3- 
2 keV) timing data. For the latter, the availability of 
spectral information and the sufficient photon statistics 
permit us to consider two energy bands: 0.3-0.7 5 and 
0.7-2 keV, where the cool and hot thermal component 
dominate, respectively. Since the effects of the system 
geometry and strong gravity are purely achromatic, the 
availability of multiple energy bands is very useful as it 
enables a test of the validity of the atmosphere model, 
which exhibits characteristic energy dependent effects 
(as described in §2.2). Due to the 10-fold better pho- 
ton statistics of the XMM-Newton dataset, compared to 
the ROSAT PSPC and Chandra HRC-S data, it pro- 
vides the best constraints on the desired parameters. We 
have considered the other data, nonetheless, as they pro- 
vide a useful consistency check of the model. The best- 
fit parameters and confidence intervals were determined 
by manually searching the x 2 hyperspace for the mini- 
mum. In principle, an additional strong constraint can 
be imposed by requiring that the derived temperatures 
and radii from these fits be within 3er of those obtained 
from the spectral analysis. Unfortunately, due to the 
limited photon statistics as well as the relatively large 
cross-calibration uncertainties of the detectors, which are 
clearly evident in Figure 4, the 3er range of derived tem- 
peratures (Ti = 1.35 - 1.95 MK and T 2 = 0.39 - 0.60 
MK) is not particularly constraining. 

The atmosphere model is able to reproduce the shape 
of the pulse profile and variation in pulsed fraction with 
energy remarkably well, as is evident in Figures 5 and 6. 
As a consequence of the numerous free parameters and 
the complex covariances that exist between them, ac- 
ceptable values of x 2 were obtained for a relatively wide 
range of parameter values. Fortunately, the strong corre- 
lations can be weakened if we fix the angle ( to its most 
likely value. Radio timing observations indicate that the 
orbital i nclination of the J0437-4715 binary syste m is 
i = 42° ([van Straten et alJl200lt iHotan et al.ll2006h . It 
is highly probable that the spin axis of the MSP is closely 
aligned with the orbital anguar momentum vector as a 

5 For this energy range, the 0.44—0.51 keV band was excluded to 
eliminate an instrumental feature specific to the fast timing mode 
of the EPIC-pn detector. 
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Fig. 5.— (Top) XMM-Newton EPIC-pn 0.3-2 keV pulse profile 
of PSR J0437-4715. The solid line shows the best fit curve for 
an off-center dipole field and a H atmosphere with R = 9.6 km, 
a = 30°, C = 42°, Aa = -12°, A</> = -20°, Ti = 2.02 MK, 
Rl = 0.2 km, T 2 = 0.48 MK, R 2 = 3.5 km, and N H = 2 X 10 19 
cm -2 . The dotted line is of a blackbody that best fits the spectrum 
of J0437— 4715 for the same assumed compactness and geometry. 
The dashed line corresponds to the best fit H atmosphere model for 
a centered dipole. (Bottom) Best fits to the XMM-Newton EPIC- 
pn pulse profiles in the 0.3—0.7 keV and 0.7-2 keV bands (upper 
and lower curve, respectively) assuming an off-center dipole. The 
choice of phase is arbitrary. 

direct result of the spin-up process during the LMXB 
phase. If this is indeed the case, the angle between the 
spin axis and the line of sight must be C ~ 42°. This 
constraint significantly constricts the allowed parame- 
ter space and thus improves the accuracy with which 
the other parameters can be determined. In particu- 
lar, for the fit to the XMM-Newton EPIC-pn pulse pro- 
files we find the ±1ct ranges to be: R — 6.9 — 10.6 km, 
a = 25° - 90°, Aa = -50° - 20°, A0 = -(23° - 14°), 
Ti = 1.4-1.85 MK, T 2 = 0.4- 0.54 MK, R 1 = 0.1-0.36 
km, and i?2 = 2.0 — 3.5 km. All ranges quoted repre- 
sent lcr confidence intervals. Varying iVn over the range 
(1 — 3) x 10 19 cm" 2 does not cause significant changes 
in the derived results since the interstellar extinction to- 
wards this MSP is neglegible above ^0.3 keV. The de- 
rived values of A(f> and Aa translate into an offset of the 
magnetic axis from the center of the star of Air = 0.8 — 3 
km (lcr). For R we are also able to derive 90% and 
99.9% limits of 6.8 — 13.9 km and >6.7 km, respectively 
(see Fig. 7); similar limits are obtained even without 
the constraint on (. The fits to the ROSAT PSPC and 
Chandra HRC-S pulse profiles are generally consistent 
with these results if we account for the systematic uncer- 
tainties encountered in the spectral fits. We note in pass- 
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Fig. 6.— Chandra HRC-S (top), and ROSAT PSPC (bottom) 
pulse profiles of J0437-4715. The solid lines show the best fit 
model to the XMM-Newton EPIC-pn data. The flux normalization 
R^ff/D 2 has been adjusted to correct for the discrepancies in the 
absolute calibration of the different detectors (see text for details). 



ing that the best fit parameters for PSR J0437-4715 de- 
rived bv lZavlin fc Pavlov! (fl998h using the ROSAT PSPC 
as well as their assumption of a centered dipole do not 
yield acceptable fits to the Chandra HRC-S and XMM- 
Newton EPIC-pn pulse profiles, even if we account for 
cross-calibration uncertainties. 

We have found that a and Aa are the two parameters 
that are most strongly covariant with R (see Fig. 8). 
Aa, in particular, closely mimics the effects of M/R on 
the observed pulse profiles (Fig. 3), thus weakening the 
constraint on R. If a and Aa for this MSP can be deter- 
mined by independent means (e.g., from radio polariza- 
tion measurements), the accuracy of the measurement of 
R can be greatly improved. 

We have conducted the entire fitting procedure out- 
lined above for a BB emission model as well. We find 
that the thermal portion of the X-ray spectrum is well 
described by a two-temperature BB model. This is ex- 
pected, given the close qualitative similarities between 
the BB and H atmosphere spectra. However, a BB 
model is incapable of reproducing the pulse shape and 
pulsed fraction for any R. Furthermore, a BB cannot 
account for the observed changes in the pulsed fraction 
as a function of photon en ergy (see Fig. 5) . Ther efore, 
we confirm the findin g of iPavlov &: Zavlinl (|1997f ) and 
IZavlin fc Pavlov! (|1998f ) that a blackbody is not a valid de- 
scription of the surface of 30^31-^115, implying that an 
optically thick, light-element atmosphere must be present 
on the stellar surface. 
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Fig. 7. — The effect of the stellar radius on the x 2 statistic of the 
fit to the XMM-Newton EPIC-pn pulse profile of PSR J0437-4715 
for M = 1.4 Mq. The solid line corresponds to a fixed value of 
f = 42°, while the dotted line is for no constraint on f. The dashed 
lines show the 68%, 90%, 99%, and 99.9% confidence levels for one 
interesting parameter. 
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Fig. 8. — 68%, 90%, and 99.9% confidence contour plots showing 
the correlation between the stellar radius R and a (top) and Aa 
(bottom) for M = 1.4 M in the fit to the XMM-Newton EPIC-pn 
pulse profile of PSR J0437-4715. 



4. DISCUSSION 

The principal motivation of our analysis is to extract 
information regarding the compactness of the NS, mag- 
netic field configuration, and radiative properties of the 
stellar surface. It is important to realize, however, that 
the method of inferring these properties is inherently 
model dependent. Thus, it is necessary to examine the 
validity of the approximations and assumptions we have 
made as well as any potential bias that they may intro- 



duce in the derived results. This is especially critical for 
determination of M/R if this approach is to be used to 
place stringent constraints on the NS equation of state 
(EOS). 

One of the underlying assumptions of the atmosphere 
model is that the source of heat is below the atmospheric 
layer. At first, this seems to be at odds with the predic- 
tion of pulsar electrodynamic models that the MSP polar 
caps are heated by a return current of ultra-relativistic 
electrons/positrons from the magnetosphere above the 
polar caps. However, the penetration depth of the im- 
pinging particles is significantly larger than the charac- 
teristic depth of the hydrogen atmosphere. In particu- 
lar, the stopping column of a relativistic electron in hy- 
drogen, assuming energy loss via bremsstrah lung and a 
negligible magnetic field, is y ~ 60 g cm~ 2 (Tsai 1974, 
and references therein). For comparison, in the case of a 
hydrogen atmosphere with logT c ff ~ 5.6 — 6. 4 the charac- 
terist ic depth is ~0.01-1 g cm -2 (cf Fig. 2 of lZavlin et alj 
1996). Therefore, the bulk of the energy of the returning 
particles is deposited well below the depth at which the 
atmosphere becomes optically thick. 

For PSR J0437-4715, the relatively long spin period 
(P = 5.76 ms) suggest that the rapid rotation has a neg- 
ligible effect on the shape of the star. Hence, our assump- 
tion of spherical Schwarzschild geometry is appropriate. 
Also, based on the pulsar period (P — 5.76 ms) and 
acceleration-corrected period derivative (1.86 x 10 20 s 
s _1 ) we find that the magnetic field strength at the sur- 
face of the MSP polar caps is B surl cx (PP) 1 / 2 w 3.3x 10 8 
G so th e assumption of an unmagnetized atmosphere is 
correct (|Zavlin et al.lll995a[) . 

In our analysis of PSR J0437-4715 we have ignored 
any thermal emission outside of the polar cap regions. 
Hubble Space Telescope FUV observations of this MSP 
(Kargaltsev et al. 2004) indicate that the rest of the star 
is at a temperature T e g ~ 10 5 K. Given that this emission 
peaks in the UV, its contribution to the X-ray flux is 
insignificant. 

Another simplifying assumption we have introduced is 
that of point-like emission regions. To determine the 
effect of this approximation, we have compared model 
pulse profiles of point-like and extended circular emis- 
sion areas with R c g ~ 2 — 4 km, comparable to those 
inferred for the cooler emission component of the spec- 
trum of PSR J0437-4715. We find that the maximum 
difference between the two cases, occurring at (j> w and 
4> w 0.5, is only ~1%, which is significantly lower than 
the statistical uncertainty of the presently available tim- 
ing data. Therefore, the hot spots are truly point-like 
so this approximation is more than adequate. Note that 
for X-ray timing data with better photon statistics, the 
exact size and shape of the emission regions may become 
important and may, in principle, be measurable quanti- 
ties. 

We have also assumed that the two thermal compo- 
nents occupy the same (0, a) location on the NS sur- 
face. The lack of any appreciable phase shift between 
the cool (0.3-0.7 keV) and hot (0.7-2 keV) components, 
as determined from the highest quality (XMM-Newton 
EPIC-pn) data, implies that the two emission regions are 
effectively concentric. The detection of any such phase 
difference would provide interesting information regard- 
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ing the temperature distribution across the polar caps. 

It is important to emphasize that the possible presence 
of faint PL emission below ~2 keV in the spectrum of 
PSR J0437-4715 does not prohibit reliable constraints 
on the properties of the underlying NS. This compo- 
nent only introduces a small increase in the uncertainty 
of our measurement. With future deep phase-resolved 
spectroscopic observations the true nature of the PL can 
be determined and the temporal properties of this com- 
ponent can be properly modeled, resulting in improved 
constraints on the desired NS properties. 

Along with insufficient photon statistics, one of the 
main factors limiting the accuracy of our measurement 
of M/R is the uncertainty in the calibration of the X- 
ray detectors used in this analysis. As noted in §3, this 
problem precludes a precise determination of T\ and T 2 
from the spectral fits, which can provide much tighter 
constraints on M/R than obtained in our analysis. Note 
that the uncertainty in detector performance has lesser 
impact on the fits to the pulse profiles, which are inte- 
grated over a relatively wide energy band so that devia- 
tions in the spectral response tend to average out. More- 
over, this uncertainty mostly effects the flux normaliza- 
tion R 2 e / D 2 of the pulse profile and has little bearing on 
the validity of our measurement of R and the conclusions 
regarding the presence of an atmosphere on the NS sur- 
face and magnetic field geometry. Improvement in the 
detector calibration will enable much tighter constraints 
on M/R. 

5. CONCLUSIONS 

We have presented a model of thermal emission from 
rotating NSs covered with an unmagnetized hydrogen 
atmosphere, applicable to systems for which the radi- 
ation is localized in small regions on the surface. Such 
a configuration is expected in MSPs. Indeed, the model 
lightcurves (Figs. 1 and 3) are qualitatively similar to the 
observed thermal pu lse profiles of nearby MSPs in the 
field of the Galaxy (|Becker fe Achenbachl l2002t IZavlinl 
2006). In addition, the broad modulations and moder- 
ate pulsed fractions (<50%) of the lightcurves produced 
by our model are in full agreement wit h the limits ob- 
taine d for the thermal MSPs in 47 Tuc (| Cameron et al.l 
12001 . 

An application to the nearest radio MSP J0437-47I5 
shows that the spectral and temporal properties of its X- 
ray emission are easily reproduced by our model, without 
the need to resort to contrived assumptions regarding 
the properties of t he NS atmosphere (composition and 
thickness, see, e.g., IHo et al.ll2007llMori fc HoH 20071 and 
the temperature distribution across the stellar surface. 
This excellent agreement allows valuable insight into im- 
portant NS properties that are not measurable by other 
means. Specifically, we confirm that the radiative prop- 
erties of the NS surface are fully consistent with a light 
element atmosphere model and are poorly reproduced by 
a blackbody. This provides compelling evidence that the 
surface of this MSP is indeed covered by a gaseous layer 
of hydrogen. Note that this result does not necessarily 
apply to all varieties of NS systems given that MSPs are 
formed via a special evolutionary channel, involving pro- 
longed accretion of gas from a close stellar companion. 
The pulse profile of J0437-47I5 suggests that the global 
magnetic field configuration closely resembles a dipole, 



albeit one significantly offset from the stellar center. It 
is unclear whether this displacement is primordial, i.e., 
the result of the core-collapse that formed the NS, the 
result of magnetic field burial during the LMXB phase 
or a consequence of magnetic fie l d evol ution and migra- 
tion as proposed by iRudermanl (|1991[ L We note that 
small scale deviations from a pure dipole, with charac- 
teristic scales <CR, may still exist near the NS surface. 
Ultra-deep X-ray spectroscopic and timing observations 
may, in principle, reveal such "microstructure" , which 
could be manifested in the form of temperature varia- 
tions across the face of the polar caps. 

Perhaps the most important practical application of 
our model is a constraint on the compactness of the NS, 
which has important implications for determination of 
the true NS EOS. In particular, the restriction of R > 6.7 
km (at 99.9% confidence) for a M — 1.4 M© star im- 
plies that PSR J0437-4715 is not ultra-compact, i.e. not 
smaller than its photon sphere. For this and other bi- 
nary MSPs a very accurate, independent measurement 
of the NS mass is possible using radio timing observa- 
tions. Thus, in principle, a unique determination of the 
stellar radius can be achieved. Note that for J0437- 
4715, the two mass mea surements of M = 1.58 ± 0.18 
(|van Straten et al.ll200l"l) and M = 1.3 ± 0.2 M Q (Hotan 
et al. 2006) currently available do not provide very tight 
constraints on the allowed masses. 

Despite the inherent model dependence of this ap- 
proach and the relatively large number of free param- 
eters we are able to extract valuable information regard- 
ing the poorly understood properties of pulsars and neu- 
tron stars, in general. Morover, as discussed in §4, the 
conclusions drawn from our analysis arc fairly robust 
in the sense that they are weakly sensitive to the un- 
derlying assumptions and approximations. Future deep 
observations of PSR J0437-4715 and other MSPs can 
provide even more stringent constraints on fundamental 
NS paramete rs. One close MSP (D k , 300 pc), PSR 
J0030 +0451 ([Becker fc Achenbachl [2001 lLommen et al.l 
2006), exhibits a double peaked X-ray pulse profile, indi- 
cating a viewing angle and pulsar obliquity substantially 
different from those of J0437-4715. As such, this pulsar 
is ideal for an independent verification of this method. 
An upcoming deep XMM-Newton observation will allow 
a detailed study of PSR J0030+0451. 

MSPs are much better suited for this analysis than 
other NS systems such as X-ray binaries, isolated NSs, 
and normal pulsars. In the latter objects there are nu- 
merous complications arising due to the effects of the 
strong magnetic field on the emergent thermal radia- 
tion, the unknown temperature distribution across the 
surface, severe reprocessing of the thermal radiation by 
the magnetosphere, and the uncertain em ission altitude 
abov e the NS surface (e.g. in X-ray bursts. ICottam et alJ 
2006). Their low magnetic fields, point-like emission 
regions, and non- variable, dominant thermal emission 
make MSPs suitable laboratories for tests of fundamental 
NS physics and constraints on their EOS. 
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